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^ ■ Abstract 

a: .. 

We determine the renormalization group invariant quark mass corresponding to the 
sum of the strange and the average light quark mass in the quenched approximation of 
QCD, using as essential input the mass of the K-mesons. In the continuum limit we 
find (M s + M)/i*K = 0.874(29), which includes systematic errors. Translating this non- 
perturbative result into the running quark masses in the MS-scheme at /J, = 2 GeV and 
using the quark mass ratios from chiral perturbation theory, we obtain m s (2GeV) = 
97(4) MeV. With the help of recent results by the CP-PACS Collaboration, we estimate 
that a 10% higher value would be obtained if one replaced Fk by the nucleon mass to 
set the scale. This is a typical ambiguity in the quenched approximation. 

June 1999 



^PARC Advanced Fellow 



1 Introduction 

Quark masses are fundamental parameters of the standard model, which have to be 
determined from experimental observations confronted with theoretical predictions [jl|]. 
At present the most precise theoretical predictions which allow for the determination 
of ratios of the three light quark masses are based on chiral perturbation theory || . A 
detailed analysis yielded || 

M u /M d = 0.55 ± 0.04 , M B /M = 24.4 ± 1.5 (1.1) 

with 

M = i(M u + M d ). (1.2) 

Unlike these ratios, the overall magnitude of the quark masses is not accessible to chiral 
perturbation theory combined with experimental data alone. One should therefore 
determine a particular linear combination of quark masses by comparing lattice QCD 
predictions J4|-|l6|] or QCD sum rules [17-27] to experiments. 



In this work we use the masses of the K-mesons and a computation in the quenched 
approximation to QCD to determine M s + M. Our analysis employs the 0(a) improved 
lattice theory, the quark mass is renormalized completely non-perturbatively [gg] and 
the continuum limit is taken (with a rate proportional to a 2 ). Hence, with respect to 
the last two points, it improves on many previous calculations. 

In general quark masses are scale- and scheme dependent quantities. It is therefore 
desirable to compute the renormalization group invariant quark masses, which - being 
both scale- and scheme-independent - are naturally taken as fundamental parameters 
of QCD. We recall that they are defined in terms of the high energy behaviour of the 
running masses fn(fi): 

Mi = lim {(2&o5 2 (/i)r do/26o m 4 Gu)) , (1.3) 

6 = ll/(4vr) 2 , do = 8/(4vr) 2 . (1.4) 

On the other hand, the renormalization group invariant masses Mj are related to the 
bare current quark masses mi by 

Mi = Z M rm (1.5) 

with a (flavour independent) renormalization factor Zm which was recently computed 
by the ALPHA Collaboration [p8| , |29| . This non-perturbative result is the basis of our 
present calculation] 2 ] The renormalization problem and its solution were discussed in 

2 Since the bare mass is involved, Zm depends on the regularization. The complete calculation of |28| 
was done in O(a) improved lattice QCD, which we use here as well. 



detail in references [28,29]. The current quark masses themselves, are defined through 
the PCAC relation, 

df.A^x) = (mi + mj)P(x), (1.6) 

in terms of the axial current, 

A^x) =A{x)l l j.l5'ipj{x), (1-7) 

and the pseudoscalar density, 

P(x)=^(x) 75 ^(x). (1.8) 

Applied to the vacuum-to-K matrix elements it reads 

M s + M = Z M (m s + m)=Z M ^ L ml, (1.9) 

where -Fk is the K-meson decay constant, and Gk denotes the vacuum-to-K matrix 
element of the pseudoscalar density.P] Eq. (03) is the fundamental relation which we 
shall exploit in this work. It is used in the following way: we compute Zyi-^ and 
multiply with the experimental squared mass, 



m 



K = i(m| + +t4o)qcd = (495 MeV) 2 , (1.10) 



to obtain M s + M. By the subscript "QCD" we indicate that we have used the masses 
in pure QCD with electromagnetic interactions switched off, since obviously the lattice 
QCD result is valid for a world where a em = 0. In practice this is achieved by sub- 
tracting an estimate of the electromagnetic effects from the experimental numbers. The 
numerical estimate in eq. ( |1.10| ) was obtained from Dashen's theorem [30] being well 



aware that the accuracy of this estimate may be only around 0.5% in mj. [ plj| . 

The combination Z^-^r- carries the dimension of an inverse mass. Therefore it is 
necessary to choose another dimensionful observable to form a dimensionless ratio which 
has a continuum limit (this is often called "setting the scale"). Choosing r$ |32| or Fk 
for this second observable and extrapolating to the continuum limit yields the results 
quoted in the abstract. 

In the following we shall first explain our strategy to deal with some technical 
difficulties in the computation of Fk/Gk- Since it involves a computation of the ratio 
.Fpg/Gps for mass-degenerate mesons, we discuss in Sect. || the dependence of various 
observables on the difference of quark masses m-j — rrij. In Sect. || we show that finite 
size effects are negligibly small in our calculation. The main result for the quark mass 
is presented in Sect. ||, where, among other issues, the extrapolations to the continuum 
are discussed. We then proceed to estimate the ambiguity which originates from the 

3 Our convention is that Fk, Gk denote the matrix elements of the bare operators. Renormalization 
factors are written explicitly. For the O(a) improved lattice theory they can be found in the appendix. 



fact that the quenched approximation can not be expected to describe the real world 
properly. As a byproduct, we present the calculation of the kaon decay constant and 
masses in the vector channel in the continuum limit. We finish with a discussion of 
our results and some comments on how the method may be extended to determine 
coefficients of the chiral Lagrangian. 

2 Strategy 

2.1 Chiral perturbation theory 

Let us first recall what chiral perturbation theory can predict concerning the quark 
masses [0,Blp3,34]. Chiral perturbation theory is based on nothing but the very general 



assumption that chiral symmetry is broken spontaneously in the limit M = M s — ► 0. 
This allows for a quantitative description of the pseudoscalar sector in terms of a low- 
energy effective Lagrangian, where quark masses appear as "kinematical variables" just 
like the energy in scattering processes. The parameters in the chiral Lagrangian are 
independent of the quark masses. At order p 4 , most of the parameters have been deter- 



mined by comparison to experimental data (see e.g. [35] for numerical values). However 



the coefficient of the chiral symmetry breaking quark mass term in the Lagrangian (de- 



noted by B in |34|), cannot be determined from experimental data alone. It can only be 
fixed when one particular quark mass is known. Since all parameters are independent 
of the quark masses, this may be done at a convenient reference point. It is important 
to realize that this reference point does not have to correspond to a physical quark 
mass. The procedure can then be extended to determine other parameters in the chiral 
Lagrangian more precisely, once additional observables are known for suitable quark 
masses and with sufficient accuracy. All of this can potentially be achieved by lattice 



QCD calculations and can improve predictions such as eq. (LI). Also the justification 
for the truncation of chiral perturbation theory can and should, of course, be checked. 
A particularly convenient way of applying this idea in practice is as follows. We 
introduce the ratio 

Fps 
R(mi,mj) = -r- , (2.1) 

such that 

mi + rrij = R(rrii, rrij) mf>g(mj, rrij) . (2-2) 

Defining a reference quark mass 

m P s(m re f , m re f) = m K , (2.3) 

with the kaon mass already discussed in the introduction, the ratio R for the physical 
quark masses may be written as 

R(m,i, rrij) = T(x i ,x j )R(m veU m re f) , x t = mi/m vei (2.4) 



with a function T(xi,Xj) which can be computed in chiral perturbation theory. In 
particular, using (3^, [3^] one finds 

T(m s /m. re f,m/m r ef) ~ 1 and thus 2m ro f ~ m s + m . (2-5) 

The corrections to the above equation are small, but the overall uncertainties associated 
with this statement require a detailed investigation of the phenomenology. The reason 
is that the errors of the parameters in the chiral Lagrangian are correlated, which makes 
it non-trivial to estimate uncertainties. Nevertheless we expect eq. ( |2.5| ) to be correct 
to within about 10%. 

In Sect. H we shall follow a complementary approach, by investigating pseudoscalar 
meson observables as a function of the average quark mass and the difference of the 
quark masses. We find that the dependence on the latter variable is rather small. It 



then follows trivially from our definition of T that eq. (2.5) is valid rather precisely. 

Before turning to that investigation let us briefly explain why our strategy is ad- 
vantageous. The reasons for quenched QCD and full QCD are somewhat different but 
related, having both to do with the difficulties of extrapolations to the region of small 
quark masses which is at present inaccessible to direct Monte Carlo calculations. 

2.2 Quenched approximation and full QCD 

The lattice calculations presented in this paper have been performed in the 0(a) im- 



proved theory 36] with non-perturbative improvement coefficients |$7| and mass-renor- 
malization factor, Zyi p3]. Details pertaining to, for instance, 0(a) correction terms in 
the currents can be found in the appendix. 

To motivate our strategy let us first review the straightforward approach for com- 
puting light quark masses. First one assumes isospin symmetry of the ratios R, i.e. 
effects of 0(m u — m^) are neglected, which is well justified given the smallness of this 
mass difference. Then one chooses an overall scale, say the hadronic radius 7*0, and 
determines the bare light quark mass and strange quark mass such that r^m^ and 
rom-K agree with the experimental numbers. Except for the bare coupling, all parame- 
ters in the QCD Lagrangian are then fixed and the current quark masses are given by 
eqs. (|2l|,|2T2D, with 



(0|^(0)|K(p)) = i Pli F K , (0|P(0)|K(p)) = G K . (2.6) 

Here |K(p)) denotes a pseudoscalar state with momentum p, standard infinite volume 
normalization and the proper flavour quantum numbers. After renormalization the 
quark masses can be extrapolated to the continuum. 

While this approach will ultimately be a clean and straightforward way to determine 
the quark masses, it poses two problems in the quenched approximation. 

• The quenched approximation is expected to be misleading for very light quark 
masses such as rh. A particular indication of the failure of the quenched approxi- 
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Figure 1: The functions Q(m,y). Qm denotes the case of pseudoscalar masses and Qf their decay 
constants. The average of rrii and rrij is between w 1.4m IG [ and « 2.4m re f. 



mation in this regime is the presence of logarithmic terms in chiral perturbation 
theory which have no counterpart in the full theory 



With our fermion action and lattice spacings, we cannot perform calculations for 
quark masses which are below approximately half of the strange quark mass. For 
smaller masses, the Dirac operator has unphysical zero modes on a certain fraction 
of configurations contributing to the path integral ("exceptional configuration") 







The second problem is of a more technical nature. It can in principle be circumvented by 
choosing a suitable action, but then also Zm has to be recomputed. However, the first 
point must be taken seriously, since we want to obtain results which are not misleading 
with respect to the real theory of interest, full QCD. 

It is furthermore advantageous to use the same strategy in full QCD. The reason 
is that it is very likely that for some time to come, lattice simulations will only reach 
down to quark masses somewhere around m = m ro f(~ m s /2). The most precise way to 
extrapolate further is then given by chiral perturbation theory as discussed above. 

3 Non— degenerate quarks 

In this section we investigate quantitatively how the pseudoscalar mass and decay con- 
stant depend on the difference of the two quark masses in quenched QCD. Apart from 



some consistency checks which were performed with the data of the ALPHA Collab- 
oration, the numerical results presented here were obtained during the course of the 
UKQCD simulations [Q using the 0(a) improved action. 

We consider the observables of interest as a function of the average current quark 
mass m = — ^ — 3 - and the parameter y defined by 

y = { Xi - Xj) 2 = (rrii - mj) 2 /m 2 ei , (3.1) 

with TO re f as determined later in Sect. [|. In particular, we want to show that the 
functions 

Hu(m, y) = to ps and H F (m, y) = F PS [1 + b A am q ] (3.2) 

have little dependence on y. In other words, the ratios 

Qu(rn,y) = Hu(m,y)/Hu(m,0) , (3.3) 

Q F (m,y) = H F (m,y)/H F (m,0), (3.4) 

are close to unity. Note that the two functions Qm and T are equivalent at the special 
point to = TO re f , since QM^ref > y) = 1/T(1 + \ y/y, 1 — \ y/y)- In order to investigate the 



y-dependence numerically, we start from the results for amps and a-Fps of UKQCD [42] 
for three values of the quark masses at each of the two bare couplings (3 = 6.0, 6.2. Three 
combinations with rrij ^ rrii are then available for each of the two corresponding lattice 
spacings, 0.09 fm and 0.07 fm. The only numerical difficulty in examining the ratios 
Qx is that the denominators Hx(jn, 0) are in general not available directly. They may, 
however, be replaced by an interpolation 

H M (m,0) = hfm + hfm 2 + hfm 3 , (3.5) 

H F (m,0) = hi + h\ to + h\ to 2 , (3.6) 

with h^ determined from the observables computed for the three available degenerate 
mesons (toj = rrij). One has to check that these interpolations are stable. This was 
tested by taking different subsets of three mass points of the data in Table [l], extracting 
hjf and comparing the resulting functions Hx{m, 0). Deviations between these different 
interpolations were found to be negligible compared to the statistical precision of the 
ratios Q themselves. 

Before discussing the results, let us mention one technical point in the numerical 
analysis. For convenience, the above procedure was carried out using the subtracted, 
improved bare quark masses, m q = m q (l + b m am q ) , m q = too — to c , instead of the 
current quark masses themselves .0 For the question addressed in this section, this 
makes no difference because to = Z m fh q holds up to the usual 0(a 2 ) errors, with some 
mass- independent renormalization factor Z m . 

4 Flavour indices are suppressed, here, and we use 6 m = —1/2 — 0.0962go- 



p 


K 


am 


amps 


amy 


aFps 


aGps 


6.0 


0.1335 


0.0466(1) 


0.3884(10) 


0.5289(26) 


0.0969(4) 


0.5997(31) 




0.1338 


0.03856(6) 


0.3529(11) 


0.5077(31) 


0.0943(4) 


0.6022(34) 




0.1340 


0.03311(7) 


0.3275(11) 


0.4938(36) 


0.0924(4) 


0.6015(39) 




0.1342 


0.02759(8) 


0.3001(12) 


0.4804(44) 


0.0905(4) 


0.5978(45) 


6.1 


0.1342 


0.04161(3) 


0.3306(5) 


0.4542(17) 


0.0870(3) 


0.7454(22) 




0.1345 


0.03305(3) 


0.2947(5) 


0.4328(23) 


0.0841(3) 


0.7476(26) 




0.1347 


0.02734(3) 


0.2687(6) 


0.4188(28) 


0.0820(3) 


0.7458(30) 




0.1349 


0.02160(4) 


0.2399(7) 


0.4056(39) 


0.0799(3) 


0.7405(37) 


6.2 


0.1347 


0.03241(3) 


0.2683(5) 


0.3748(21) 


0.0743(3) 


0.8887(33) 




0.1349 


0.02661(3) 


0.2430(6) 


0.3588(26) 


0.0721(3) 


0.8908(38) 




0.13515 


0.01934(3) 


0.2080(6) 


0.3388(39) 


0.0693(3) 


0.8857(48) 




0.1352 


0.01788(3) 


0.2004(6) 


0.3348(43) 


0.0687(3) 


0.8830(51) 


6.45 


0.13485 


0.02477(2) 


0.1975(8) 


0.2749(25) 


0.0533(5) 


1.261(10) 




0.1351 


0.01734(2) 


0.1650(9) 


0.2558(37) 


0.0503(4) 


1.265(14) 




0.1352 


0.01439(2) 


0.1505(9) 


0.2480(46) 


0.0490(4) 


1.261(16) 




0.1353 


0.01142(2) 


0.1347(10) 


0.2398(59) 


0.0477(4) 


1.250(19) 



Table 1: Results for masses and unrenormalized (ratios of) matrix elements in all simulation points 

In Fig. [l] we show our numerical results]^] They show that Qm has a dependence on 
y which is below the level of a percent, and Qp is seen to be independent of y to within 
our statistical precision of around 0.3%. It should be kept in mind that Fig. |] tests the 
y-dependence for m/m re f = 1.4-2.4, while later we shall be interested in kaon physics, 
where m/m re { ~ 1. It remains a (plausible!) hypothesis that Qx are close to unity 
also for this lower value of m. Nevertheless, our investigation quantifies the smallness 
of y-effects for the first time and suggests that the correction terms in eq. ( [2.5| ) are 
negligible. Based on this analysis we shall not distinguish between 2m re f and m s + m 
from now on. 

4 Finite size effects 



Our numerical results, which are listed in Table ffl, have been obtained for approximately 
constant physical volume, T x L 3 , where T = 2L and L ~ 3ro = 1.5fm. We used exactly 



the numerical methods described in [43], and in particular systematic errors due to 



5 The numerical values for (m q ) re f were taken from the analysis described in Sect. 



excited state contributions were checked to be small compared to the statistical errors. 
Further details about the simulation are described in the appendix. Since the precision of 
the results is quite good one has to investigate whether they may be affected by the finite 
size of the system at the level of their statistical accuracy. The most relevant observable 
for the determination of the quark masses is the ratio R = (mi + Trij)/mpo. Since rrii 
is defined through the PCAC relation, it is independent of the volume |36|] (apart from 
small lattice artefacts of order a 2 ). We are therefore predominantly interested in the 
volume dependence of the pseudoscalar masses but will also consider the decay constant 

Fps- 

In the pseudoscalar sector the leading finite size effects can be reliably calculated 
in chiral perturbation theory. The reason is that for vanishing quark masses the pseu- 
doscalar mesons are Goldstone bosons. Their interactions become weak for small en- 
ergies and these are responsible for the leading finite size effects for large but finite 
volumes. 

Gasser and Leutwyler have reported results for the two-point function of the axial 
current at space-like separations [44]. This result, given for finite temperature and in 



finite volume is easily adapted to the situation which is of interest here, namely the 
pseudoscalar mass on an L 3 -torus (temperature zero) mps(L), and the corresponding 
decay constant Fps(L). The resulting formulae are 

HS^L i = l^ 9W+ 0(e-^),, = mpsi (4.1) 

mps(oo) N{ F£ s 

m^r-1 = -N t ?k g (z) + 0(e-&), (4.2) 

-Fps(oo) F£ s 

9{z) = A-, r%- z2x - 1/(4x) = A- Ki(z), (4.3) 



7T 2 Z 2 Jo x<1 27I" 2 , 

where K\(z) denotes a modified Bessel function. Here a comment is in order. This is a 
result of the first non-trivial order in chiral perturbation theory in full QCD. A priori 
one cannot expect it to be accurate if the pseudoscalar meson masses are too large. On 
the other hand, one knows that the finite size effects discussed here are of order e~ z , as 



long as the pseudoscalar is the lightest particle in the theory [ 45 , 46 1 . It is plausible that 
the prefactor of the exponential is not dramatically different for heavier pseudoscalar 
mesons. We may therefore take the above equations as a reasonable estimate of finite 
size effects. Being interested in the order of magnitude of these effects, we also do not 
worry about the difference of the quenched approximation and full QCD. 

In our numerical calculations (see Tables [I] and |3]), we are in the range z = mpsL > 
4.3 and consequently eqs. fl4.lD -(fI~2J) predict corrections which are below the level of 
0.5% for Fps and 0.1% for mps- In order to check this estimate, we have done some cal- 
culations on lattices which are larger than our standard L ~ 1.5 fm. Entirely consistent 
with the above formulae, we found no significant changes in Fps,mps, for instance 



r O m PS 


P 


7 R 


M B +M 
(Fps)k 


n>(-Fps)R 


myr 


(F PS )n 

my 


1.5736 


6.0 


0.1939(30) 


0.810(11) 


0.3746(45) 


2.428(32) 


0.1543(26) 




6.1 


0.2077(28) 


0.824(10) 


0.3952(47) 


2.457(32) 


0.1609(26) 




6.2 


0.2160(30) 


0.843(11) 


0.4020(48) 


2.361(40) 


0.1701(32) 




6.45 


0.2205(46) 


0.851(19) 


0.4070(60) 


2.441(71) 


0.1667(52) 




CL 


0.2300(69) 


0.874(27) 


0.4146(94) 


2.386(66) 


0.1761(75) 


3.0 


6.0 


0.1966(26) 


1.450(17) 


0.4063(47) 


2.637(22) 


0.1541(20) 




6.1 


0.2104(26) 


1.471(17) 


0.4288(49) 


2.666(20) 


0.1609(20) 




6.2 


0.2181(28) 


1.501(18) 


0.4363(51) 


2.606(25) 


0.1675(23) 




6.45 


0.2236(37) 


1.518(28) 


0.4421(62) 


2.677(42) 


0.1652(32) 




CL 


0.2324(57) 


1.553(42) 


0.4507(99) 


2.645(41) 


0.1709(49) 



Table 2: Extra-/interpolations in rQm PS and resulting continuum limits (CL) including all errors 



13 = 6.2, amps = 0.208 : 



mps(24a) 
mps(32a) 
Fp S (24a) 
Fps(32o) 



- 1 = 0.003(4) (ss 0.0009) , (4.4) 

1 = -0.003(8) (« -0.003) , (4.5) 



where the numbers in parentheses are the estimates from the above equations with 
JVf = 2. 

We conclude that even with the numerical precision of Table |l| finite size effects 
are negligible for the quark masses considered. Of course, as described by the above 
formulae, finite size effects grow rapidly when the pseudoscalar mass becomes smaller 
and larger volumes would be necessary for quark masses which are significantly smaller 
than the ones we used. 



5 Degenerate quarks, continuum limit 

We can now proceed to compute the desired quantity tq(M s + M). The first step is to 
evaluate R for degenerate quarks as a function of ?~o m PS- Relegating all details of the 
exact numerical procedure, such as the definition of m q and the different improvement 
coefficients, to the appendix we directly show the mass dependence of R in Fig. ^a. It 
is apparent that R is almost constant as a function of the quark mass (or r^mpg). It 
is therefore easy to extrapolate to the desired point ?"o m PS = r o m K = 1-5736, slightly 
outside of the range where we have numerical results. The extrapolation is performed 
linearly in ?~o m PS> usm § the three closest data points. Simply taking for instance the 
closest data point or an average of all of them would change the final result for the 
quark mass by a negligible amount. 
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Figure 2: Mass depend ence and extrapolations for the two smallest values of the lattice spacing. The 
ratio R is defined in eq. (A..6). 



The second step is now to form the combination 



r {M s + M) = Z M - 



R\ 



'crops' 



=1.5736 



x 1.5736 



m 



(5.1) 



collecting all errors, including the ones on tq [47] and the /3-dependent one of Zm [28]. 
We then extrapolate to the continuum limit linearly in a 2 . As seen in Fig. ^a, the 
dependence on the lattice spacing is significant - in contrast to the quark mass on 



smaller volume [48]. As a safeguard against higher order terms in the lattice spacing, we 
therefore include only the three points with smallest lattice spacing in the extrapolation 
and obtain our main result 



r (M s + M) = 0.362(12). 



(5.2) 



10 



This result contains the /3-independent part of the uncertainty in Zm of 1.3% [28]. 
At this point one may be concerned about the validity of an a 2 -extrapolation, since 
6a — bp is mostly known only in perturbation theory. However, this combination is 



found to be tiny at the one-loop level |53| and turns out to be rather small also non- 
perturbatively |50|,51|. We have checked explicitly that such small values of 6 a — 6p do 
not affect our continuum extrapolation. 

In the whole procedure one may replace R, defined in eq. ( |2.2| ), by the ratio of 
the current quark mass and the pseudoscalar mass squared. One then ends up with 
slightly different lattice spacing effects and somewhat different statistical errors. We 
have performed also that analysis, using the current quark mass averaged over a range 
of timeslices around the central one in our lattice. The results differ by up to 2% from 
the ones presented above for finite lattice spacings but are indistinguishable after the 
continuum extrapolation. 

Instead of the quark mass in units of ro, we may also compute the combination 
(M s + M)/(i ? K)R, where (-Fk)r is the kaon decay constant. As its experimental value we 
take (-Fk)r = 160(2) MeV. We note that electromagnetic effects cannot be subtracted 
in analogy to eq. ( 1.1CD , chiefly because -Fko is not known experimentally. We emphasize 



that this is truly an alternative procedure which amounts to computing the combination 

M^ + M M 1 

x 1.5736, (5.3) 



(-Fk)r m Z P rfiGps[l + b P am q 



-j<mjg=1.5736 



where the ratio M/m and the renormalization constant Zp have to be taken at a common 
scale for which we resort to the value fj, = l/(1.436ro) chosen in [28|. Discretization 



effects are clearly expected to be different in this case: they are known to be non- 
negligible in the product ro(i ? ps)R p%- A non-perturbative estimate being unavailable 



at present, we have to rely on the perturbative value 6p = 1 + 0.153 <?o [ 53 1 . Extrapo- 
lations in r^mps and the continuum extrapolations are shown in Fig. ||b and Fig. yb. 
Concerning the former, the mass-dependence of l/{Gps[l + 6pam q ]} is stronger than 
the mass dependence of R, but within our errors it is perfectly linear and the extrapola- 
tion is easily done. In comparison to the quark mass in units of ro, we find discretization 
errors which are roughly only half as big. 

In order to confirm that it is legitimate to perform an extrapolation linear in a 2 , 
even though the improvement coefficient 6p is known only perturbatively, we have re- 
peated the whole analysis after setting 6p to its tree- level value 6p = 1. Results after 
extrapolation changed by much less than a percent. Our final continuum result is 

^=0.874(29). (5.4) 

For illustration we may translate to physical units, setting ro = 0.5 fm [47| and (-Fk)r = 
160(2)MeV. We thus obtain 

M s + M = 143(5) MeV from eq. (5.2) , (5.5) 

M s + M = 140(5) MeV from eq. (§3|) . (5.6) 
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Figure 3: Continuum limit extrapolations of several observables. Full symbols show the extrapolated 
values. Dashed lines represent the extrapolation function, which are continued outside the fit range 
towards larger lattice spacings. 

As will be discussed in more detail below, this assignment of physical units is ambiguous 
in the quenched approximation. One should be well aware that the solid results are given 
in eq. ( fT^ ) and eq. ( |5.4| ). 

For future reference we list various dimensionless quantities both for ?~o m PS = 
1.5736 and for r^mpg = 3 in Table |2[ Our results before and after the continuum 
extrapolation can be found in that table. 

6 Ambiguities in the quenched approximation 

In the previous section, we have mainly used tq as our reference scale. However, in 
general one expects that by choosing different experimental inputs one will also get 
somewhat different results when working in the quenched approximation. Here we 
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would like to quantify this ambiguity for our determination of M s + M. 

A first estimate is obtained by computing the combination tq(Fk)r, using our 
results for the pseudoscalar decay constant including the non-perturbative renormal- 
ization factor Z\, and comparing it to its experimental value. Of course, on the basis 



of the results for M s + M presented in eqs. (5.5) and fl5.6|) one may not expect any 
significant deviation. However, the following detailed analysis may still be instructive, 
in particular since it yields ro(i 7 k)R hi the quenched approximation. 

The procedure to extract tq(Fk)r is entirely analogous to what was done in the 
previous section. The extrapolations in r^mpg (for mass-degenerate quarks) and the 
continuum extrapolations are shown in Fig. [|a and Fig. |3|c, respectively. In the contin- 
uum limit we find 

quenched: {r Z A F PS [l + b A am q ]} r 2 m 2 =L5736 = 0.415(9) (6.1) 

compared to 

experiment: 0.5 fm x (F k )r = 0.405(5) (6.2) 

with the experimental value of (-Fk)r- Clearly the result obtained in the quenched 
approximation agrees well with experiment when the approximate value ro = 0.5 fm is 
employed. However, it would be premature to conclude that the ambiguity is small, 
since the agreement may be special to the case of (-Fk)r- 

Another alternative to set the scale is provided by the nucleon mass, m^. Compared 
with other widely used quantities like m p , the nucleon has the advantage of being a 
stable hadron. Although we have not computed win ourselves, we can still use the 
precise published results obtained by the CP-PACS Collaboration (l(| and combine 



them with the values of r^ja reported in [47|. We obtain 



quenched: toton ~ 2.6 (6-3) 

compared to 

experiment: 0.5 fm x m^ = 2.38 , (6-4) 

which represents a 10% difference. Of course also the result for M s + M in MeV would 
change by the same amount if one used win instead of tq or (-Fk)r as experimental 
input. This may serve as a rough estimate of the inherent ambiguity in the quenched 
model. 

In principle one could also use directly the results from recent comparisons of the 



quenched light hadron spectrum with experiment 16,54] to estimate this ambiguity. In 
ref. |16| statistically significant deviations of up to 11% are observed when the scale is 
set by m p , which, as a resonance, has a fairly large width. Since we consider it safer 
to set the scale by the nucleon mass, we have refrained from quoting this number as 
the typical uncertainty for the case at hand. However, using the published numbers 
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Figure 4: Mass dependence of pseudoscalar decay constant and vector meson mass. 



of [16], it is not difficult to estimate the quenched results for the other hadron masses 
when one chooses m^ as input - at least within a precision of 2-3%. Interestingly it 
then turns out that the masses of the stable hadrons agree with experiment to within 
about 4%. On the other hand, for unstable hadrons one observes differences which can 
be as large as their widths. This may not be too surprising, since resonance effects are 
not controlled in the lattice calculation. In short, it is difficult to assess the relevance 
of the deviations observed in the quenched hadron spectrum for our analysis, and thus 
we stick to our above estimate of around 10% for the ambiguity in question. 

Finally, let us discuss the quark mass dependence of the flavour non-singlet vector 
meson masses in more detail. First we note that effects of the differences of quark 
masses can be shown to be unimportant in the same way as in Sect. |3] and we restrict 



our attention to the masses of mass-degenerate mesons (rrii 



rm 



Our aim is to map 



out the quark mass dependence of my in the continuum limit: we pick certain values of 
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Figure 5: Quark mass dependence of flavour non-singlet vector meson masses. Experimental masses 
are shown as asterisks and their width is indicated by a dotted line. 



q = m/m s = M/M s , and for each of these values and for each lattice spacing we then 
perform an inter-/extrapolation of my as a function of the quark mass to determine 
my(m = qm a ). Here, m s = 1.921 m re f (see eq. (|1.1|) ) is used and m re f is known from 
the previous section. At fixed value of q we then extrapolate r^my to the continuum 
limit including our data for all values of the lattice spacing, after observing that the 
lattice spacing dependence is very weak. 

In Fig. H we plot the continuum results as a function of q. They are compared to 
rriK*, mp and m^. Since all of these states are resonances, we also indicate their width 
in the figure. Additional reservations apply to the inclusion of the <j) meson in this 
comparison. It is a mixture of flavour octet and singlet components, while the lattice 
calculation is for a pure octet state with mass-degenerate quarks: the disconnected 
quark diagrams occurring in the singlet channel are not accounted for, and mixing with 
glueballs is neglected. Although all these effects may be argued to be irrelevant in the 
quenched approximation, the experimental value in Fig. S does include them and it is 
somewhat surprising that our values for my at q = 1 come so close to m^. 

If we were to ignore these problems and determined M s from the requirement that 
r$my at that quark mass agrees with rrjm^, it is evident from Fig. || that we would 
obtain a result for M s which is compatible with the numbers quoted in the previous 
section. However, the error would be considerably larger. Further continuum results 
such as romx*, shown also in Fig. ||b and Fig. ||d, and (-Fk^r/ttik* are listed in Table ^ 
for future reference. 

To summarize, we have directly seen ambiguities of about 10% owing to the use of 
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the quenched approximation. This has to be kept in mind when quark masses in MeV 
are quoted. Still we have demonstrated that within the quenched approximation a total 
precision of 3% (see eq. ( |5.2| ) and eq. (fO])) can be achieved. 

7 Discussion 

Noting that ratios of light quark masses can at present best be calculated in chiral 
perturbation theory, we discussed a strategy to compute the overall scale of the light 
quark masses and applied it to the test-case of quenched QCD. 

The problem has essentially two parts. First, the renormalization problem should 
be solved non-perturbatively. Here we have used the recent solution of the ALPHA- 
collaboration [28|, 29 j which connects the bare current quark masses to the renormaliza- 



tion group invariant ones. Second, one particular quark mass should be computed by 
matching with pseudoscalar meson masses. It is not necessary to perform this matching 
at physical values of the quark masses, because the parameters in the effective chiral 
Lagrangian are mass-independent. For both the quenched approximation and full QCD 
it appears convenient to use two different (flavours of) quarks with equal mass M re f 
such that raps is equal to the kaon mass tok- 

The ratio of M re f to the physical quark masses can then be computed in chiral per- 
turbation theory. In Sect. || we have shown that pseudoscalar masses, tops, have very 
little dependence on the difference of the two quark masses in the quenched approxi- 
mation. From this one has 2 M re f ~ M s + M, which is also true in chiral perturbation 
theory (in full QCD!). 

By making the identification 2M re f = M s + M we have computed (M s + M)tq and 
(M s + M)/(Fk)r, incorporating all systematic errors (such as finite size effects), and 
have taken the continuum limit of these quantities. These solid results are given in 
eq. (|J) andeq. |[j ). 

A conversion of the results to physical units is necessarily ambiguous, since quenched 
QCD does not describe the real world. Depending on the choice of quantity some- 
what different results in MeV are obtained. It is interesting to note that using either 
ro = 0.5 fm or (-Fk)r = 160 MeV gives indistinguishable results at our level of precision 
of 3%. However, this does not mean that M s + M has to be very close to our quenched 
result also in full QCD. Indeed, we have estimated that replacing (-Fk)r by ton (or 
other stable light hadron masses) would give quark masses in MeV, which are larger by 
about 10%. 

Usually the running quark masses in the MS scheme, to^^), are quoted. Their 
relation to the renormalization group invariant quark masses, which we computed here, 
has been discussed in Sect. 2 of |28fl . For convenience, we include a table with the 
perturbative conversion factors m^jg(/x)/M computed by numerical integration of the 
renormalization group equations with the n-loop approximations of the renormalization 
group functions (/?, r in the notation of [p3]) for n = 2,3,4. One should be aware that 



the result Ai-^ = 238(19) MeV [28 enters the numbers of Table |3[ The uncertainty in 
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H [GeV] 




^Ms(^)/ M 






2-loop 


3-loop 


4-loop 


1.0 


0.80279 


0.83585 


0.84449 


2.0 


0.70388 


0.71830 


0.72076 


4.0 


0.64079 


0.64880 


0.64981 


8.0 


0.59549 


0.60055 


0.60105 


90.0 


0.49937 


0.50105 


0.50112 



Table 3: Factors to convert the renormalization group invariant mass into the MS scheme at scale /j, 



A translates into 1.5% in m^((j,)/M at \i = 2 GeV and 2.5% at \i = 1 GeV. If desired, 
Table || may be used to estimate frly^(fi). A typical result combining eq. (|5.6P, eq. (1.1) 
and the table is 



m s (2 GeV) = 97(4) MeV with 4-loop running in the MS scheme . 



(7.1) 



This illustrates that a high level of precision can be reached by state of the art lattice 
techniques. 

Effects of dynamical fermions may first be examined in the theory with unphys- 
ically large quark masses. We have also given results for a reference mass which is 
roughly twice the strange quark mass. One should quantify what the effects of dynam- 
ical fermions are at such a point, where a good accuracy may be achieved, and then 
move on to the more chiral regime. 

Also for some quantities not related directly to quark masses, precise results have 
been obtained. They refer to the continuum limit of the quenched approximation and are 
summarized in Table 0. Most notably, decay constants were computed with a precision 
not much worse than the experimental one. This was achieved using the method to 
compute hadronic correlation functions proposed in p3[ . In this context we should 
also comment on the continuum extrapolations in the O(a) improved theory. They are 
certainly necessary! For the quark mass in units of ro, the difference between its value 
at a lattice spacing a = 0.1 fm and at a = amounts to about 15%. Other quantities 



show smaller lattice spacing effects as was observed also in a finite volume study [48]. 
Although clearly one would have hoped for a weaker a-dependence, we note that their 
order of magnitude is not much different from what is known already for pure gauge 



theory observables [47]. 

Finally let us come back to the role of lattice QCD in determining parameters of the 
chiral Lagrangian. Figs, y and ||a show to what precision the quark mass dependence 
of observables in the pseudoscalar sector can be computed. Such information, once 
available in full QCD, will allow to reduce the errors in the chiral Lagrangian. For 
instance, we have checked that just the mass independence of R puts strong restrictions 



17 



on these parameters. As usual, estimates of the uncertainties is a delicate issue and a 
detailed investigation of the potential of this approach will be left for future work. 



The project of computing light quark masses was started quite a while ago by 
the ALPHA collaboration. Here we presented the final results. The basis of our ap- 
proach was developed together with M. Liischer, S. Sint and P. Weisz. We would like 
to thank them for a most enjoyable collaboration. In particular, R.S. and H.W. thank 
M. Liischer for numerous enlightening discussions. Furthermore, we appreciate useful 
correspondence with G. Ecker and D. Wyler on chiral perturbation theory and electro- 
magnetic corrections. 

The numerical computations were performed on the APEIOO/Quadrics computers of 
DESY Zeuthen and on the Cray T3D at EPCC Edinburgh. We thank these institutions 
for their support. This work was in part supported by EPSRC grant GR/K41663 and 
PPARC grants GR/K55745 and GR/L29927. 

A Numerical details 

A.l Improvement coefficients, renormalization factors 

Our calculations have been performed using the 0(a) improved Wilson action defined 
in ref. [37], which can be consulted for any unexplained notation. In particular, the 



improvement coefficients c sw and c\ were taken from eqs. (5.15) and (6.4) of that refer- 
ence. 

The renormalized axial current and pseudoscalar density for quark flavours i and j 
are defined as 

(A R ) / ,(x) = Z A (l + 6 A am q )(Ai) M (x), (A.l) 

P R (x) = Zp(l + b P am q )P(x), (A.2) 

where P{x) = tp i (x)'j5^j(x), and the improved, unrenormalized axial current is given 
by 

(Ai) M (s) = ^(xh^s^Or) + ac A ± {d; + fy) P(x). (A.3) 



The renormalization factor Z\ was calculated in J49J , and here we have used its param- 
eterization in eq. (6.11) of that paper. 

The renormalization factor Zm, which relates the current quark mass in the 0(a) 
improved theory to the renormalization group invariant quark mass, has recently been 



determined [28|. In that paper also the scale- and scheme-dependent factor Zp has been 
computed in the Schrodinger functional (SF) scheme at a fixed scale of L = 1.436 tq. 
Here we use their representation in terms of polynomial fit functions, viz. 

Z M (go) = 1-752 + 0.321 (J3 - 6) - 0.220 {(3 - 6) 2 , (A.4) 
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Z P (go,L/a) L= iA36r = 0.5233 - 0.0362 ((3 - 6) + 0.0430 (/? - 6) 2 , (A.5) 
(3 = 6/gl 6.0 < (3 < 6.5. 



As explained in Subsect. 6.2 of ref. [28|, the uncertainty in Zu is split into a /3-dependent 
part of 1.1%, which enters any continuum extrapolation, and a /3-independent error of 
1.3% which must be added (in quadrature) to the extrapolated result. The typical 
accuracy of Zp as given by the above polynomial is 0.5%. 

The factor Zyi renormalizes the lattice counterpart of the ratio R(mi,mj) defined 
in eq. (|2.1|). In the O(a) improved theory it is given by 

R(mi,mj) = — — [1 + (6a - &p)am q ] , (A.6) 



where 



I {K) 4 + K),} = 1.(1 + 1. 1). (A.7) 



and Kj, Kj are the hopping parameters of flavours i and j. 



The combination 6a — 6p has been found to be small in perturbation theory [53] 



and also non-perturbatively [|50|, 51|. A method how to determine the individual coef- 
ficients 6a and 6p non-perturbatively has been proposed in ref. [55|, but no results for 
our choice of action have been reported so far. 

The values for the critical hopping parameter k c have been taken from Table 1 of 
ref. |J7|]. When necessary they have been interpolated linearly to the desired (3- value. 

A. 2 Hadronic correlation functions, reference scale 

Our calculation of hadronic correlation functions using Schrddinger functional boundary 
conditions follows exactly the procedures outlined in a previous paper E3]. In particular, 
this reference contains the definitions of the correlation functions fp, / A , ky and /i, 
which we have computed to extract hadron masses and matrix elements for pseudoscalar 
and vector mesons. 

Estimates for meson masses were obtained by averaging effective masses, 

am eS , x (x + f ) = In (f(x )/f(x + a)) , X = PS, V, (A.8) 

over a suitably chosen interval £ m i n < xo < t max . Similarly, the ratio -Fps/Gps was 
obtained by averaging the combination 

Gps mps Jp(x ) 
and the individual matrix elements Fp$ and Gpg are extracted from 



-Fps - 


= 2(m PS L 3 r 1 /2 e (zo-T/2) mps WOJ ; 

v7i 


(A.10) 


Gps = 


= 2(m PS /L 3 ) 1 /V a;o - T / 2 ) mp s /Pjfo) ; 


(A.ll) 
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where again an average over the interval [t m irn^max] is taken. The averaging intervals 
used here were the same as those in Table 1 of ref . [03] . 

We have also calculated an unrenormalized current quark mass defined by 

. v U d o + 9o)/a(xo) + c A ad%d fp(x ) 

m(x ) = — . (A.12) 

2/p(x ) 

The values for the current quark mass listed in Table [l] have been obtained by averaging 
m(xo) over the interval T/4 < xq < 3T/4. 

For the hadronic radius rQ which was used to set the scale, we have taken its 
parameterization quoted in [E7J, viz. 

ln(a/r ) = -1.6805 - 1.7139 ((3-6) + 0.8155 {(3 - 6) 2 - 0.6667 ((3 - 6) 3 , (A.13) 

accounting for the uncertainty as quoted in [47]. In addition to the results obtained at 



(3 = 6.0 and 6.2, which have been published in [43|, we have also computed correlation 
functions at [3 = 6.1 and 6.45. A brief summary of our simulation parameters including 
-^rneas, the number of gauge field configurations for which the correlation functions were 
computed, is presented in Table |4|. 



(3 


T/a 


L/a 


L/r 


1 'mcas 


6.0 


32 


16 


2.981(12) 


1000 


6.1 


40 


24 


3.795(17) 


800 


6.2 


48 


24 


3.261(16) 


800 


6.45 


64 


32 


3.060(17) 
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Table 4: Simulation parameters 
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